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£N| ■ Abstract. This paper deals with an Gierer-Meinhardt model, with three sub- 

stances, formed Reaction-Diffusion system with fractional reaction. To prove 
global existence for solutions of this system presents difficulties at the bound- 
ednees of fractionar term. The purpose of this paper is to prove the existence 

^sj of a global solution using a boundary functionel. Our technique is based on 



(N 

< 



- '-I 



the construction of Lyapunov functionel. 



1. Introduction 



In recent years, systems of Reaction-Diffusion equations have received a great 
deal of attention, motivated by their widespread occurrence in models of chemical 
and biological phenomena. These systems are divided into celebrated classes; sys- 
tems with dissipation of mass and systems of Gierer-Meinhardt. In this paper we 
deal with this last. 

£> ' In the study of the various topics from plant developmental; Meinhardt, Koch 

and Bernasconi [10] proposed Activator-Inhibitor models (an example is given in 
section 5) to describe a theory of biological pattern formation in plants (Phyl- 
^1 , lotaxis). 

• ' We assume a Reaction-diffusion system with three components: 

O; f fjf - ai A« = /(tt, v, to) = a - ftitt + vn ff, +c) 

(1.1) { ^-a 2 Av = g(u,v,w) = -b 2 v + J^ *eft,t>0, 

ff - a 3 Aw = h(u, V, w) = -b 3 w + ^g^ 



with Neummann boundary conditions 
du dv dw 



d-2) ararw =0 on ^ x{f>0} ' 

and the initial data 

(u(0,x) = ip^x) > 
v(p,x)=<p 2 (x)>0 onfi, 
w(0 1 x) = ip 3 (x) >0 

and tpi € C (ft) for all i = 1, 2, 3. 

Here ft is an open bounded domain of class C 1 in R , with boundary 9ft and 

d 

— — denotes the outward normal derivative on 9ft. 
or] 
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c, ,pi, qi, r,: are non negative with a, 6j, a% > 0, indexes for all i = 1, 2, 3. 

(1.4) <pi - 1 < max^p 2 min ( — , — , 1 ) , p 3 min — , — , 1 

I V.92 + 1 r 2 J V r 3 + ! 93 

Put Aij = °' aj for all i,j — 1,2, 3. Let a, (3 and 7 be positive constants such 



M 3 — 2 x /a-a 

that where 



a > 2max<j 1, 2 ^ 3 } , (1.5) 



61 

-> CM?„ 

/3 



> 2^? 2 , (1.6) 



and 



2 



(L7 > {w- A ")(^- A °>)>{^ A ™- A " A « 

The main result of the paper reads as follows: 

Theorem 1. Suppose that the functions f,g and h are satisfing condition (1-4)- 
Let (u (£, .) , v (t, .) , w {t, .)) be a solution of (l.l)-(l.S) and let: 

(1.8) L(t) = f "° ( *' X 1 ~J x. 

Then the functional L is uniformly bounded on the interval [0,T*],T* < T max . 
Where T max (Huolloo > Ikolloo > 1 1 i^o 1 1 00 ) denotes the eventual blow-up time. 

Corollary 1. Under the assumptions of theorem 1 all solutions of (l.l)-(l.S) with 
positive initial data in C (U) are global. If in addition b\, 62, b%, a > 0, then 
{u, v, w) are uniformly bounded in £1 x [0, 00) . 

2. Previous Results 

The usual norms in spaces L P (J1), L°°(f2) and C(fl) are denoted respectively by: 

IK = ±-Jju(x)\ p dx; (2.1a) 

IML = max|u(aO|, (2.1b) 

IMIc(H) = max|w(a;)|. (2.1c) 

In 1972, following an ingenious idea of Turing. A [17], Gierer.A and Meinhardt. H 
10^ proposed a mathematical model for pattern formations of spatial tissue struc- 
tures of hydra in morphogenesis, a biological phenomenon discovered by A.Trembley 
in 1744 [TBJ.lt is a system of reaction-diffusion equations of the form: 

(2.2) (f-«xA« = cr-MU+S for all s e fl ,i >0 

I H - a 2 Av = -vv + ^ 

with Neummann boundary conditions 

on r)n 

(2.3) — = 0, and — = 0, x € dtl, t > 0, 
on or) 
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and initial conditions 

(2.4) (^'S^fwn' * efi > 

where f2 C K is a bounded domain with smooth boundary 9f2, oi, 02 > 0, /i, z^, o > 
0,p, g, r and s are non negative indexes with p > 1. 

Global existence of solutions in (0, 00) is proved by Rothe in 1984 [14] with 
special cases: p = 2, g = 1, r = 2, s = and iV = 3. The Rothe's method cannot 
be applied (at least directly) to the general p, q, r, s. It is desirable to consider 
the p, q, r, s originally proposed by Gierer-Meinhardt. Wu and Li [8] obtained the 
same results for the problem (2.1)-(2.3) so long as u, w _1 and a are suitably small. 
Mingde, Shaohua and Yuchun llj show that the solutions of this problem are 
bounded all the time for each pair of initial values if 

r s + 1 

^— ^ < 1. (2.5b) 

r 

Masuda. K and Takahashi. K 9 we consider a more general system for (u, v) : 

( oa\ / F ~ aiAu = CTl ^ ~ /iM + Pl &' U ^ W 

[ -> \ &L-a 2 Av = a 2 (x)-vv + p 2 (x,u)^ 

with 0-1,0-2 6 c 1 (o) , 0-1 > 0,0-2 > o,pi,p 2 g c 1 fn x s+^) n i°° fn x 1+) 

satisfying pj > 0, p 2 > and p, q, r ,s are nonnegative constants satisfying (2.5a). 
Obviously, (2.4) system is a special case of (2.6) system. In 1987, Masuda. K and 
Takahashi. K [5] extended the result to 2=- < j^ under the condition a\ > O.In 
2006 Jiang.H [7] under the conditions (2.5a) - (2.5b) , (pj_,(p 2 £ W 2 ' 1 {&) ,1 > 
max{iV, 2} , ^ = ^ = on d£l and ip t > 0, <p 2 > in IT.Thcn (2.6) system has 
a unique nonnegative global solution (u,v) satisfying (2.3)-(2.4). 

3. Preliminary Observations 

It is well-known that to prove global existence of solutions to (1.1) — (1.3) (see 
Henry [6]), it suffices to derive a uniform estimate of ||/(u, v, w)\\ , \\g (u, v, w)\\ p 
and \\h(u, V, w)\\ on [0;T max [ in the space L p (f2) for some p > iV/2.0ur aim is 
to construct polynomial Lyapunov functionals allowing us to obtain L p — bounds 
on u; v and w that lead to global existence. Since the functions /, g and h are 
continuously differentiable on M^_, then for any initial data in C(O), it is easy to 
check directly their Lipschitz continuity on bounded subsets of the domain of a 
fractional power of the operator 

aiA 

(3.1) O = - I a 2 A 

a 3 A 

Under these assumptions, the following local existence result is well known (see 
Friedman [3] and Pazy [T3]). 

Proposition 1. The system (1.1) — (1.3) admits a unique, classical solution (u; v; w) 
on (0,T max [xfi. 

(3.2) IfT max < 00 then lim (||u (t, .)\\ 00 + \\v (t, OIL + \\w (t, ,)\\J = 00. 

* / J- m ax 
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Remark 1. This proposition seems to be well-known (Dan Henry 6 ). Neverthless 
we could not find it in the literature in the form stated here and stated in the book 
of Franz Rothe ( |14j pp: 111-118 with proof). Usually the explosion property (3.2) 
is only stated for some norm involving smoothness, but not the Loo — norm. 

4. Proofs 

For the proof of theorem 1, we need a preparatory Lemmas, which are proved in 
the appendix. 

Lemma 1. Assume that p, q, r, s, m, and n satisfying 

p — 1 / q m 

< mm — — , — , 1 

r \s + 1 n 

For all h, I, a, j3, 7 > 0, there exist C — C (h, I, a, (3, 7) > and 9 = 9 (a) G (0, 1), 
such that 

T p~l+a T r+a / a \ e 

(4.1) a—, — — <(3 ,,,. - + C — — ) , x>Q,y>h,z>l 

Lemma 2. Let fi,T > and fj = fj(t) be a non-negative integrable function on 
[0, T)and < 9j < 1 (j = 1, ..., J). Let W = W(t)be a positive function on [0, T) 
satisfying the differential inequality 

(4.2) ^M < -»W (t) + £/,- (*)***' (*) , 0<t<T. 

Then, we obtain that 

(4.3) W(t)<K, 0<t<T, 

where k is the maximal root of the following algebraic equation: 

(4.4) x - V ( sup / e-"(*-«/i(0de > ) ^ =W(0). 

-—{ \0<t<TJo J 



Moreover, ifT= +00, then 



limsupM^ (t) < k c 

t/'oo 



where Koo is the maximal root of the following algebraic equation: 

x-V] (limsup f e-rt-Qfjffldt] X 6 * =0. 

t— J y t/'oo Jo J 

Lemma 3. Let (u (£, .) , v (t, .) , w (t, .)) be a solution of (l.l)-(l.S), then for any 
(t,x) in (0,T max [xfl we get 

u(t,x) > e~ blt mm(tp 1 (xj) > 0, 
(4.5) { v(t,x)>e- b2t min(ip 2 (x))>0, 

w(t,x) > e~ bst min(cp 3 (x)) > 0. 
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proof of theorem 1. Differentiating L (t) with respect to t yields 



L'(t) = 



'■^ d * u - P^+T^ d ^ ~ ^^TT dtW dx ' 



replacing dtu, dtv and dtw with its values in (1.1), we get 

r ( 
L'(t) 



a ia ^Au - a 2 /3^^Av - a^^^Aw 



,,P2+° 



\ I n y ' " a v.* 

\ ' v*l +Pw1 (lU r l +C) ^ 1 ,92+/3+l u ,''2+T 

/+ J, 



"f V q 3 +P t 



„J>3 + <* 



da;, 



r-3+T+l 



era 



'v^uF 



where / contains laplacian terms and J contains the other terms 

,a—\ 



I 
J 



a\a 



U: 



V W 1 



Audx — a 2 j3 



ir 



yP+T-Wl 



(-b ia + b 2 /3 + b 3 -f)L(t) 

f u pi+«-i 
+a / -= — =-. rdx 

,a-l 



Avdx — a 3 j 



u P2+a 



If 



n 



—Awdx, 



(!»' 



<72+^+l ?/ ,r 2 +7 ' /„ w ?3+/3 u) r 3 +7+l 



II 



P3+a 



dx 



-act 



a vPw"f 



dx. 



Starting with estimation of /: 

Using Green's formula for terms f Q ^^Audx , f n v ^ lwl Avdx and J n „^° +1 Awdx 
we get 

-a x a (a - 1) ^ |Vu| 2 + a t af3 jl'^ VuVv + maj J^ VuVw 
+a 2 paJ^VuVv - a 2 /3 (/? + 1) ^^ \\7v\ 2 - a,^ v , + f w ^ +1 VvVw^ 
+a 3 ja v f w ^ +i \7u\7w - a 3 j/3 v0+ ^ +i VvVw - a 3 j (7 + 1) ^° +2 \Vw\ 2 




v 



a-2 



/3+2 w7 +2 



{QT) ■ T 



dx, 



where 



-a/3 ^i±^ 



/3 7 ^±^ 



aia (a — 1) 

Q - ! -a/3^p a 2 /3 (0 

_ a7 £ilM ^ 7 £2±£3 037 ( 7 +l) 

<5 is matrix of quadratic form compared to: vwVu, uw\7v and uuVw, which is 

written in the matric form:T = (uwVw, uw\7v, uv\7w) . 

Q is definite positive if, and only if all its principal successive determinants are 

positive. To see this, we have: 

1. Ai = 01a (a - 1) > 0. Using (1.5), we get Ai > 0. 

^ Z{JVl) =« 2 Ai a2 (^^-^ 2 ). Using (1.5) 

and (1.6), we get A 2 > 0. 

3. Using theorem 1 in S.abdelmalek and S. Kouachi [T] we get (a — 1) A3 = 

(a - 1) \Q\ = a (a 7 /3) 2 ai a 2 a 3 ((^^ - A\ 2 ) (*£*&- A{ 3 )-(^A 23 A 12 A 

Using (1.5)-(1.7), we get A 3 > 0. 



2. A, 



dx. 



1:? 



) 2 )- 
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Consequently we have / < 0, V (£, x) e [0, T*\ x ft. 

Now we estimate J : 

J = (-ha + b 2 p + 637) L (t) 

r u pi+a— 1 r u P2+a f U P3+a 

-i-Q; / dx — (3 1 dx — t / dx 

J Q Vl 1+ Pw'l(w ri +C) ' J a v4 2+ P +1 W r2 +"l ' Jq t;93+/3y,r- 3 +7+l 

f u °'~ 1 -, 
+aa / —3 da;. 

According to the maximum principle, there exists Co dependant on (p 1 , (p 2 and ^3 3 
such that v, w > Co > 0, then we have 

q-l (3 T q-1 /3 + -f 

u «-i / u « \ — /!\- /1\~ / u Q \~ / 1 ' 

then 

^- < C 2 (-£-] " where C 2 = ( i- 

we have 

J = {-ha + b 2 p + 637) L (t) 

f u pi+a ~ 1 f u P2+a f u P3+a 
-\-(% I dx — (3 1 dx — ^y I dx 

J n V <li+(3 W ~t ( W ri 4. c ) ' ' J n u 92+/3+l w r 2 + 7 ' ^ w 93+/3 ^3+7+1 

+£ra / —5 — da;. 



Using lemma 1, V(£,x) G [0,T*] [xfi we get 



u Pi+a-l u P2+a / U« \ 8 



(4-6) " — TTi 7 ^<a — z« — 1~ < P — XTTTi ^ + C 



w <?l+/3 w 7( w ri + c ) - W 9l+/3 w 7+fl - M w 92+/3+l w r 2 +7 \vPw~t) 



or 



u pi+a-l uP3+a / U a \ 
( 4 - 7 ) « Zfl 37" < 7 Z7T ZTT + C ^ 

Using (4.6) or (4.7) then 

J < (-ha + b 2 p + b 3 i)L (t)+ c(-^—\ dx + aa C 2 (-^—^ ) d.v. 
Applying Holder's inequality, for all t in [0,T*] we obtain 

C [ ^-\dx<(((4^\ dx \((c^dx^ e 



in v^w 7 / \Jn\ vl3w ~ t J / Wf! 

then 

/ C ( "XT ) dx - ^"W- where C 3 = C M 

We have 



then 



we get 
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/ C 2 ( -Jf-7 ) " dx < C 4 L^ (t) where C 4 = C 2 |fi| 



J < (-61Q + 6 2 /3 + b 3l ) £ (*) + C 3 L e (t) + aaC 4 L— (t) , 
which implies 

J < (-ha + b 2 (3 + b 3l ) L (t) + C 5 (l° (i) + aaL^r (*)) . 

Thus under conditions (1.5), (1.6) and (1.7), we obtain 

L'(t) < (-61 a + b 2 /3 + 6 37 ) L (t) + C 5 (l 6 (t) + acrL^ (*)) , 

since — b^a + b 2 (3 + 637 < and Using lemma 2 we deduce that L(t) is bounded on 
(0,T max [ ie L(t) < k, where k dependent on ip x , ip 2 and ip 3 . □ 

Proof of Corollary 1. Since L(i) is bounded on (0, T max [ and the functions vn ^ ujri , ■, 

■^^ and ^^ are in I°°((0, T max ), L m (n)) for all m > f , then as a con- 
sequence of the arguments in Henry. D }6 j or Haraux. A and Kirane. M [3] we 
conclude the solution of the system (l-l)-(l-7) is global and uniformly bounded on 

ftx(0,+oo). □ 

5. Example 

In this section we will examine a particular activator-inhibitor model in order 
to illusrate the applicability of corollary 1 and proposition l.We assume that all 
reactions take place in a bounded domain $7 with a smooth boundary dfl. 

Example 1. The model proposed by Meinhardt, Koch and Bernasconi [10] to de- 
scribe a theory of biological pattern formation in plants (Phyllotaxis) , where u, v 
and w are the concentrations of three substances; called activator (u) and inhibitors 
(v and w) is: 

du d 2 u u .. 1 a 2 1 

(5.1) { §_ a2 |!| = _fe 2U + w 2 ; forallxen, t>0. 



dt 



M a 3 ^ = -b 3 w + u, 



Proposition 2. Solutions of (5.1) with boundary conditions (1.2) and nonnegative 
uniformly bounded initial data (1.3) exist globally. 

Proof. This model is a special case of our general model (1.1), where p\ =2,q\ = 
l,ri = 1,|>2 = 2,^2 = 0, r 2 = 0,^3 = 1,(73 = 0,r 3 = 0. These indexes realize the 
conditions of global existence: Pl ~ 1 < min ( -2W , — , 1 ) . D 

Remark 2. TTie system described by equations (5.1) exhibits all the essential fea- 
tures of phyllotaxis. 
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6. Appendix 

The purpose of this appendix is to prove lemma 1, lemma 2 and lemma 3 in 
section 4 which we have used in the proof of theorem 1. 

Proof of Lemma 1. For all x > 0, y > h, z > I we have from the inequality (4.1) 



(6.1) a- </3—^ + C, . 

and we can write 

X V ~ X P-1 / X T \ r (s + l)(p-l) n(p-l) 

a = an r a — —. — ) v q z *■ 

y q z m f ^ y s+l z nJ » 

For each e realize: < e < min ( -fr, — , 1 1 — £— ^ 



-1 / r \ -E^+e 

X F . p-i I „ X \ r I . X' \ (£±iKp-1) „ "(p-i) 



yq z m ^ \ ' yS+1 z n I \ 1t s ~^^ Z n 



yq z m \ ' n ,s+l^n I \ ,,s+l 

Then also 






(!)'(?) 



where 



S C l( ^) %+ '(^') , (6.2) 



p_l ^^i„p-i; _ pr. ( „ )( j) 



Using Young's inequality for (6.2) with taking C — C 1 r ' and 0=1 — 

where e is sufficiently small, we get inequality (6.1). □ 



a l-E=i-£ 



Proof of Lemma 2. This lemma is proved in [ Masuda.K and Takahashi. K [S], 
Lemma 2.2]. D 

Proof of Lemma 3. Immediate from the maximum principle. □ 
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